Quantum dynamics of an optomechanical system in the presence of Kerr-down 

conversion nonlinearity 
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We study theoretically nonlinear effects arising from the presence of a Kerr-down conversion 
nonlinear crystal inside an optomechanical cavity. For this system we investigate the influences 
of the two nonlinearities, i.e., the Kerr nonlinearity and the parametric gain, on the dynamics of 
the oscillating mirror, the intensity and the squeezing spectra of the transmitted field, and the 
steady-state mirror-field entanglement. We show that in comparison with a bare optomechanical 
cavity, the combination of the cavity energy shift due to the Kerr nonlinearity and increase in the 
intracavity photon number due to the gain medium can increase the normal mode splitting in the 
displacement spectrum of the oscillating mirror and reduce its effective temperature. Our work 
demonstrates that both the Kerr nonlinearity and down conversion process increase the degree of 
squeezing in the transmitted field. Moreover, we find that in the system under consideration the 
degree of entanglement between the mechanical and optical modes decreases considerably because 
of the intracavity photon number reduction in the presence of the Kerr medium. 
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I. INTRODUCTION 

In recent years, there has been an increasing interest in 
cavity optomechanical systems for a wide range of both 
experimental and theoretical investigations|lH8(.The im- 
portance of optomechanical systems is due to their po- 
tential applications in different topics of physics. They are 
promising candidates for studying quantum effects in the 
mesoscopic and macroscopic scales [9|, detection and in- 
tcrfcrometry of gravitational waves [10l4l2| and measure- 
ment of small displacements [Hj]. In a cavity optomechan- 
ical system the radiation pressure exerted by the electro- 
magnetic field induces a coupling between the intensity 
of the cavity field and the mechanical motion of a mov- 
able mirror. It is pointed out that the necessary condi- 
tion for observing quantum phenomena for the oscillat- 
ing mirror, as a mesoscopic or macroscopic object, is the 
preparation of the mechanical oscillator at low phonon 
occupancy . Although the ground state cooling of the 
mirror has not yet been achieved experimentally, it has 
been shown theoretically 0, EH EU that the ground state 
cooling of the mirror is possible in the resolved sideband 
regime where the cavity bandwith is less than mechani- 
cal oscillation frequency of the mirror. Another quantum 
phenomenon expected to be observed in optomechani- 
cal sytems is the mixing between the fluctuations of the 
cavity field around the steady state and the mechanical 
mode which leads to the normal mode splitting (NMS) in 
the displacement spectrum of the mirror [l6j|. It is well 
known from cavity QED that NMS is an absolute ev- 



' sareh.shahidani@gmail.com 



idence of strong coupling between the subsytems with 
energy exchange taking place on a time scale faster than 
the decoherence of each mode. The appearance of NMS 
in optomechanical systems is known as a direct conse- 
quence of ground state cooling of the mirror flfjf. 

One of the most important characteristics of optome- 
chanical systems is their nonlinear optical peroperties. 
The strong interaction between the cavity filed and the 
mechanical oscillations makes the cavity behave like a 
nonlinear medium since the length of the cavity depends 
upon the intensity of the field in analogous way to the 
optical length of a nonlinear material^- It has been 
shown[T3] that within the Born-Openheimer approxima- 
tion the back action of the oscillating mirror in one 
edge of an optical Fabry-Perot cavity,driven by a nearly- 
resonant laser field, induces an intracavity Kerr-like non- 
linearity. Besides this intrinsic nonlinearity, the presence 
of an optical parametric amplifier (OPA) or the op- 
tical Kerr medium [l9j inside the cavity has opened up 
a new domain for combining nonlinear optics and op- 
tomechanics towards the enhancement of quantum ef- 
fects. It has been demonstrated 18| that the presence 
of an OPA in the cavity causes a strong coupling be- 
tween the oscillating mirror and the cavity mode result- 
ing from increasing the intracavity photon number. Thus 
the enhancement of radiation pressure-induced coupling 
not only considerably improves the cooling of the me- 
chanical oscillator but also makes the observation of the 
NMS of the movable mirror and the output field mode 
accesible[18|. On the other hand, when the optomechan- 
ical cavity contains an optical Kerr medium with x 
nonlinearity, due to the photon-photon repulsion and the 
reduction of the cavity photon fluctuations, the normal 
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mode splitting weakens and the effective temperature of 
the moving mirror increases fl9j. 

The above mentiond interesting results motivated us 
towards investigation of an optomechanical system which 
contains a nonlinear crystal consisting of a Kerr medium 
and a degenerate OPA. We will show that in the system 
under consideration the competition between the increas- 
ing of intracavity photon number due to the parametric 
gain process and the reduction of the number of photons 
due to the photon blockade mechanism arising from the 
presence of the Kerr medium leads to some new inter- 
esting effects in the dynamics of the movable mirror in 
the low photon number limit and the resolved sideband 
regime. 

The paper is structured as follows. In SccHI we de- 
scribe the model, derive the quantum Langevin equations 
of motion for the system operators and hnd their steady- 
state mean values. In Scc lIIII we linearize the quntum 
Langevin equations of motion around the steady-state 
mean values, and analyze the stability conditions of the 
system. In Sec lIVI we calculate the spectrum of small 
fluctuations in the position of the oscillating mirror , the 
effective frequency and effective damping rate of the me- 
chanical oscilltor, and analyze the influence of the non- 
linear gain and anharmonicity parameter on them. In 
Sec[V]we investigate the influences of the parametric gain 
and Kerr nonlincarity on the intensity and quadrature 
squeezing of the transmitted field. In Sec I VII we examine 
the entanglement between the optical and the mechanical 
mode. Finally, we summarize our conclusions in Scc lVIII 



II. THE PHYSICAL MODEL 



the detection bandwidth is chosen such that it includes 
only a single, isolated, mechanical resonance and mode- 
mode coupling is negligible. It should be noted that in 
the adiabatic limit, the number of photons generated by 
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FIG. 1. (Color online) Schematic picture of the setup stud- 
ied in the text. The cavity contains a Kerr-down conversion 
system which is pumped by a coupling field to produce para- 
metric oscillation and induce Kerr nonlinearity in the cavity. 

Under these conditions, the total Hamiltonian of the 
system in a frame rotating at the laser frquency w^can 
be written as 



We consider a Kerr-down conversion optomechanical 
system composed of a degenerate OPA and a nonlinear 
Kerr medium placed within a Fabry-Perot cavity formed 
by a fixed partially transmitting mirror and one movable 
perfectly reflecting mirror in equilibrium with a thermal 
bath at temperature T. The movable mirror is free to 
move along the cavity axis and is treated as a quantum 
mechanical harmonic oscillator with effective mass m, 
frequency ui m and energy decay rate 7 m = to m /Q (where 
Q is the mechanical quality factor). The cavity field is co- 
herently driven by an input laser field with frequency 
and amplitude e through the fixed mirror. Furthermore, 
the system is pumped by a coupling field to produce para- 
metric oscillation and induce the Kerr nonlinearity in the 
cavity. In our investigation, we restrict the model to the 
case of single-cavity and mechanical modes [2(| [2l|. The 
single cavity- mode assumption is justified in the adia- 
batic limit, i.e., u> m << nc/L which c is the speed of 
light in vacuum and L is the cavity length in the absence 
of the cavity field. We also assume that the induced res- 
onance frequency shift of the cavity and the nonlinear 
parameter of the Kerr medium are much smaller than 
the longitudinal-mode spacing in the cavity. Further- 
more, one can restrict to a single mechanical mode when 



H = H + H 1 (1) 

where, 

p 2 1 

H Q = h(uj Q - u L )a>a + (- h -muj^q 2 ) - hg m a^aq 

Zm Z 

+ihe{a { - a), (2a) 
Hi = ihG(e w a) 2 - e~ l9 a 2 ) + hr/a^a 2 . (2b) 

The first two terms in Hq are, respectively, the free Hamil- 
tonian of the cavity field with annihilation(creation) op- 
erator a(a^) and the movable mirror with resonance fre- 
quency uj m and effective mass m, the third term describes 
the optomechanical coupling between the cavity field and 
the mechanical oscillator due to radiation pressure force, 
and the last term in Hq describes the driving of the in- 
tracavity mode with the input laser. Also, the two terms 
in Hi describes, respectively, the coupling of the intra- 
cavity field with the OPA and the Kerr medium; G is the 
nonlinear gain of the OPA which is proportional to the 
pump power driving amplitude, 9 is the phase of the field 
driving the OPA, and r\ is the anharmonicity parameter 
proportional to the third order nonlinear susceptibility 
X^ 3 ' of the Kerr medium. The input laser field populates 
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the intracavity mode through the partially transmitting 
mirror, then the photons in the cavity will exert a radia- 
tion pressure force on the movable mirror. In a realistic 
treatment of the dynamics of the system , the cavity- 
field damping due to the photon-leakage through the in- 
complete mirror and the Brownian noise associated with 
the coupling of the oscillating mirror to its thermal bath 
should be considered. Using the input-output formalism 
of quantum optics [HI, we can consider the effects of these 
sources of noise and dissipation in the quantum Langevin 
equations of motion. For the given Hamiltonian ([1]), we 
obtain the following nonlinear equations of motion 



V 

q = — , 
m 

p = -mujf n c] + hg m a^a - -f m p + £ 
a = —i(uiQ — uil)o- + igrnQO, + e 



2irja* a 2 



(3a) 
(3b) 

(3c) 



where n is the cavity decay rate through the input mir- 
ror and din is the input vacuum noise operator that is 
charactrized by the following correlation functions (25| 

< 5a m {t)5a\ n {t') > = 5(t - t'), (4a) 

< 8a ln (t)5a m (t') > = < 5a\ n {t)5a\ n {t') >= 0. (4b) 

The Brownian noise operator £ describes the heating of 
the mirror by the thermal bath at temperature T and is 
characterized by the following correlation function [24| 



< mm >- 



2tt 



-iu>(t-t') 



t)UJ 

2k B T 



' coth( — )+l]<Ax>, 

(5) 

where fcs is the Boltzmann constant. 

We are interested in the steady-state regime and in 
small fluctuations with respect to the steady state. Thus 
we obtain the steady-state mean values of p, q and a as 



0,q s 



K-iA + 2Ge l 



A' 2 



AG 2 



(6) 



(7) 



where q s denotes the new equilibrium position of the 
movable mirror and A = luq — lol — QmQs + 2?y|a s | 2 = 
A+277|a s | 2 is the effective detuning of the cavity which in- 
cludes both the radiation pressure and the Kerr medium 
effects. It is obvious that the optical path and hence the 
cavity detuning are modified in an intensity-dependent 
way. The first modification which is a mechanical non- 
linearity, arises from the radiation pressure-induced cou- 
pling between the movable mirror and the cavity field and 
the second modification comes from the presence of the 
nonlinear Kerr medium in the optomechanical system. 
It has been shown that [2(| for a given sufficiently large 
pure x^ 3 ' nonlinearity inside an optical cavity, the effect 
of photon blockade occurs as a consequence of large phase 



shifts of the cavity detuning. It has also been predicted 
[27J that the same photon-photon interaction can occur 
in the strong coupling regime in the optomechnaical sys- 
tems. In our treatment the photon-photon interaction 
due to the radiation pressue is ignorable in comparison 
with photon-photon interaction due to the Kerr nonlin- 
earity. Since A satisfies a fifth-order equation, it can 
have five real solutions and hence the system may ex- 
hibit multistability for a certain range of parameters. 



III. DYNAMICS OF SMALL FLUCTUATIONS 

In order to investigate the dynamics of the system, we 
need to study the dynamics of small fluctuations near 
the steady state. We assume that the nonlinearity in the 
system is weak and decompose each operator in Eq.© as 
the sum of its steady-state value and a small fluctuation 
with zero mean value, 



Sa, q = q s + Sq, p = p s + 5p. 



(8) 



Inserting the above linearized forms of the system opera- 
tors into Eq.([3]), the linearized quantum Langevin equa- 
tions for the fluctuation operators take the form 



d A 

—Sa 

dt 



Sp/m, (9a) 

—muj 2 n 5q + hg m (a s 5a) + a* 5a) - ^y m Sp + £(9b) 

-(iA + n)6a - ig m a* s Sq + (2Ge 10 - 2irja 2 s )Sa 1 
—Airi\a s \ 2 Sa + V2~KSai n , (9c) 
(iA - n)6a^ + ig m a s Sq + (2Ger iB + 2ir]a* s 2 )6a 
+4ir 1 \a s \ 2 Sa ji + V2^6a\ n . (9d) 

Defining the cavity field quadratures 8x = 8a + Sa^ and 
Sy = i(Sa^ — Sa) and the input noise quadratures Sxi n — 
Sa ln +Sal n and <%„ = i(5a\ n — 8a ln ) Eq.© can be written 
in the compact matrix form 



—da' 
dt 



u = Mu(t) +n(t), 



(10) 



where u(t) = (Sq,Sp,Sx,Sy) T is the vector of fluctua- 
tions, n(t) = (0,£,^/2^Sxi n ,\i / 2fi5yi n ) T is the vector of 
noise sources and the matrix M is given by 



M 



where we have defined 






1/m 










-7m 


ha+ 


—iha- 


2ia- 





-K + T 1 


Ai + £1 


2a + 





-Ai + 5i 


— K — Tl 



(11) 



a± = 9m(a s ± a*)/2, 
Ti = 2G cos 9 -ir)(a 2 -a* 2 ), 
Aj = A + 4?7|a s | 2 
5i =2Gsm6-r](a 2 s +a* 2 ). 



(12) 
(13) 
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Here, we concentrate on the stationary properties of the 
system. For this purpose, we should consider the steady- 
state condition governed by Eq. (fTD|) . The system is stable 
only if the real part of all eigenvalues of the matrix M 
are negative, which is also the requirement of the validity 
of the linearized method. The parameter region in which 
the system is stable can be obtained from the Routh- 
Hurwitz criterion [2^|, which gives the following three in- 
dependent conditions: 



si = 2kRi + cj^(7„ 



\k 1 



s 2 = {R^ 2 m - A 



-rj\a s 



>0, 



(14a) 



(14b) 



- c.c.)} > 0, 

S3 = -[RlOjf n - R2] X W^(7m + 2k) + [2kRi 

+ uj 2 n { lm + An 2 )] > 0, (14c 
where we have defined R\ and R 2 as 

R x = Aj - AG 2 + k 2 + 1277 2 |a s | 4 
+8r?Ai|a s | 2 + AiGrj(a 2 s e~ ' 



R2 



Ml 



{2r,\a s \ 4 + iG(a 2 e-* 6 



- c.c), 

C.C.)). 



(15a) 
(15b) 



IV. DISPLACEMENT SPECTRUM OF THE 
OSCILLATING MIRROR 

Since we are interested in the spectrum of fluctuations 
in position of the oscillating mirror in the presence of 
Kerr-down conversion nonlinearity, it is convenient to 
solve the time-domain equation of motion (|10[) by Fourier 
transforming it into the frequency domain. The Fourier 
transform of the time-domain operator u(t) is 



u(lS) 



1 

2^ 



(16) 



Then solving the linearized quantum Langevin equation 
for the position fluctuations of the oscillating mirror 
yields 

dq(cu) =F 1 (Lo)aLo)+F 2 (u)Sa m (u)+F 3 (Lo)6al(Lo), (17) 
in which 



md(uj) 
mayuj) 

F 3 (w) = F 2 {-uf , (18) 
with 

d(u) = (uj 2 n -u 2 + iwym)(A» + {k + iuj) 2 ) - /?!, (19) 
2hgl 



ft = ^H[ A \a s \ 2 +2G(a 2 e-'° -c.c.)}, 
m 

p 2 = -AG 2 + A 2 - Arj 2 \a s \ 4 + AiGr){a 2 s e 



c.c. 



(20) 



The first term in Eq. ([T7|) originates from the thermal 
noise while the second and third terms arise from the 
radiation pressure force . The displacement spectrum of 
the mirror is defined by 

SJuj) = -1- I dn e -^ +n)t < Sq(oj)Sq(Q)+5q(n)5q{oj) > 
An J 

(21) 

Using the following correlation functions in the frequency 
domain 

< 5a in (u)5al n (n) >= 2ttS(uj + 0), (22) 



< £(w)f(n) >= 2Trhj m oj[l+coth(-^-)]S(u;+ni (23) 

2k B ± 

the displacement spectrum of the oscillating mirror is 
obtained as 



S q {u) = h\ X \ 2 {m lm ucoMj^—) + 2 



(hg m ) 2 n 



'2k B T' m 2 \d(uj)\ 2 
\a s \ 2 5 2 +AGRc{alc-^{n + l ^)) ^ 

(fo + k 2 - uj 2 ) 2 + An 2 u 2 h K ' 

where <5f = (n 2 + uj 2 + AG 2 + A' 2 ). 

A. Normal Mode Splitting 

In order to determine the structure of the displacement 
spectrum of the moving mirror we need to determine the 
eigenvalues of the matrix iM as the solution of Eq. (fTQ| in 
the frequency domain, or the zeros of d(ui) in the denom- 
inator of coefficients in Eq. (fT8|) . First we assume that 
the radiation pressure coupling g m is zero and find the 
eigenvalues of iM as 



2 



± 



^7, 



V, 



-in± v/32- 



(25) 



Thus in the absence of optomechanical coupling the effec- 
tive frequency of the normal modes depends on /?2 which 
is a function of nonlinear coefficiens rj ,G and intracavity 
intensity . Now we assume that the optomechanical cou- 
pling g m is not zero, but the decay rates j m and k are 
negligibly small in comparison with ui m and @2 ■ Then 
the zeros of d(w) are approximately given by 



1 



4 = 



/3 2 ) 2 +4/3! 



(26) 



Hence, the Kerr-down conversion medium can alter the 
radiation pressure contribution to the NMS in two ways. 
Firstly, both nonlinearity parameters G and rj change the 
mean intracavity photon number considerably. Secondly 
the parameter /?i is a function of the gain and anhar- 
monicity parameters. 

In FigJSJ by using the experimental parameters from 
Ref . [30j . we have plotted the mean intracavity photon 
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number a s and the normalized parameter /?i//?o versus 
the normalized gain G/k for two different values of an- 
harmonicity parameter 77, where f3 — 2hg^ n (A\a s \ 2 /m) 
is the contribution of radiation pressure to the normal 
modes of displacement spectrum in the bare cavity. The 
figure shows that while there is a reduction in the num- 
ber of intracavity photons (hence the radiation pressure) 
due to the photon-photon repulsion mechanism (in the 
presence of Kerr medium) , the contribution of radiation 
pressure coupling to the normal modes which is contained 
in /?! can be increased by choosing proper values of initial 
detuning and nonlinear coefficients G, r\ and 9 in such a 
way that /3x > (3q ■ 
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FIG. 2. (Color online) (a)The mean intracavity photon num- 
ber and (b) the normalized parameter fii/Po versus the nor- 
malized gain G/k for 77 = 0.01Hz (red solid line) and 77 = 
0.05 Hz (blue dashed line) .The parameters are uj m /2-R = 10 
MHz,L = 1mm, m = lOng, Q — 5 x 10 5 Hz, input laser power 
P = 6.9 mW at A = 1064 nm, k = 0.1u m , and = tt/2 . 



It is well known that the NMS in a bare-damped op- 
tomechanical system occurs only if <?o| a s| ^ K (where 
go = \/ E/2rriui m Qm) due to the finite width of the 
peaks [l5|, [Tg|. In Fig.© we have plotted the displace- 
ment spectrum of the movable mirror below this thresh- 
old (c/ol a s| ^ K )• As expected, the radiation pressure- 
induced coupling between the oscillating mirror and the 
fluctuations of the cavity mode does not lead to NMS. 



FIG. 3. (Color online)The displacement spectrum ver- 
sus the normalized frequency uj/uj m for bare cavity (77 = 
0, G = 0)(blue thin curve), for cavity with Kerr nonlinearity 
(77 = 0.01Hz, G = 0) (orange dashed-dotted curve), with gain 
nonlinearity (77 = 0, G = 8 x 10 6 Hz) ( green dashed curve) and 
with both nonliearity(?7 = 0.01Hz, G = 8 x 10 6 #2) (red thick 
curve). The parameters are uim/2n = 10 MHz,L = 3mm, 
m = 12ng, Q — 5 x 10 Hz, input laser power P — 9 mW at 
A = 1064 nm, « = 0.02cj m , and = tt/2 . 



The presence of the Kerr medium alone dose not change 
this situation, while in the presence of the gain medium 
solely, mode splitting occurs only when the parametric 
gain G is large enough. However, as shown in Fig.([3|. 
the combination of these two nonlinearities leads to the 
appearance of NMS in the displacement spectrum. We 
can come to this conclusion from the fact that in the 
peresence of the two nonlinearities, there is a competition 
between the two quantum processes; the nonlinear gain 
medium increases the intracavity photon number and 
splits up the pump photons into two degenerate subhar- 
monic photons , while the Kerr medium shifts the energy 
levels of the cavity proportional to the anharmonicity pa- 
rameter 77 which blocks the entrance of a large number of 
photons into the cavity. The combination of the energy 
shift resulting from the Kerr medium and increasing the 
photon number due to the gain medium leads to the ap- 
pearence of NMS. It is notable that the presence of the 
gain medium decreases the degree of photon-photon re- 
pulsion due to the Kerr nonlinearity, and hence increases 
the parameter /3± in Eq. (|26p and the radiation pressure 
effect on the oscillating mirror. In Fig.Q and Fig.© 
we have plotted the displacement spectrum of the mov- 
able mirror above the threshold of splitting (po|os s | > k) 
for different values of the gain coefficient G and differ- 
ent values of anharmonicity parameter 77, respectively. 
As expected, the presence of the gain medium inside the 
cavity strengthens the coupling between the oscillating 
mirror and the cavity mode, while the presence of the 
Kerr medium inside the cavity weakens this coupling due 
to the photon-photon repulsion and shifts the modes far 
from each other. 
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FIG. 4. (Color online)The displacement spectrum versus the 
normalized frequency u)/uj m for different values of the para- 
metric gain: G = 0(green thin curve), G = 10 J Hz(blue dashed 
curve), and G = 10 7 Hz(red thick curve). The parameters are 
r] = 0.01 Hz, A = w m , 9 = 7r/2, u) m /2n = 10 MHz,L = 1mm, 
effective mass m = 10ng, Q = 5 x 10 5 Hz, P = 30 mW, 
A = 1064 nm, and k — 0.01o; m . 
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FIG. 5. (Color online)The displacement spectrum versus 
the normalized frequency u) /to m for different values of the an- 
harmonicity: rj = 0.01Hz(red thick curve) r\ — 0.03Hz(blue 
dashed curve), and rj = 0.06Hz(green thin curve) .The pa- 
rameters are G = 10 6 Hz , A = Lj m , = 7r/2, cj m /27r = 10 
MHz,L = 1 mm, m = lOng, Q = 5 x 10 s Hz, P = 30 mW, 
A = 1064nm, and k = 0.01o; m . 



Effective frequency and effective damping rate 
of the oscillating mirror 



of the Kerr and the gain nonlinearities on the effective 
frequency w e // and the effective damping rate j e ff of 
the oscillating mirror. We can find the mechanical sus- 
ceptibility of the mirror according to the dependence of 
the Sq(ui) on the fluctuations in the total force on the 
mirror JFT(qj)f35|. 



Sq(oj) = x(w)^.Ft(w). 



(27) 



According to Eq.ffTT]), SFt(uj) consists of a radiation pres- 
sure and a Brownian motion term. Thus the mechanical 
susceptibility of the mirror (x), its effective resonance 
frequency {ui e ff) 1 and its effective damping rate (7 e //) 
are, respectively, given by 



X 1 = m (^eff -U 2 ) +iuJ"/ e ff, 



2 



- 2 ) 



|ft + (- K + iW ) 2 | 2 ' 



Jeff = m(7 m + 



2f3 lK 



\fo + (-K + iw)*\ 



(28) 



(29) 



(30) 



In FigsEJa) and [6jt>) we have plotted the normalized 
effective frequency of the oscillating mirror w e f / /w m ver- 
sus the normalized frequency u)/ui m for A = uj m and for 
different values of G [Fig. ([6k,)] and different values of r\ 
[Fig.fJBb)]- As is seen, although the gain medium does 
not change the effective frequency of the moving mirror 
considerably , the energy shift of the cavity mode due 
to the Kerr nonlinearity manifests itself in the response 
frequency of the oscillating mirror. This means that the 
effective coupling between the cavity and the mechani- 
cal modes (and therefore the cooling of the mirror) takes 
place when the effective detuning of the cavity A (but 
not A) is chosen near the oscillating mirror frequency. 
In Figs(7Ja) and[7Jb) we have plotted the normalized ef- 
fective damping rate of the oscillating mirror j e f f /m"f m 
versus the normalized frequency Lo/uj m for A = u> m and 
for different values of G [Fig[7|Ja)] and different values 
of i] [FigE^b)] . As is seen, by increasing the gain pa- 
rameter the effective damping rate of the mirror motion 
increases while increasing the anharmonicity parameter 
not only leads to the frequency shift but also decreases 
the effective damping rate. 



The radiation-pressure force can modify the dynam- 
ics of the mechanical oscillator and therefore the sus- 
ceptibility of the oscillating mirror can be considered as 
the susceptibility of an oscillator , with an effective fre- 
quency and an effective damping. These quantities con- 
tain important information about the quantum behav- 
ior of the system. The radiation presure-induced change 
in the frequency of the mechanical mode is known as 
the optical spring effect [32| , and the radiation presure- 
induced change in the damping rate of the mirror shows 
how the optical field acts effectively as a viscous fluid that 
damps the mirror oscillation and cools its center-of mass 
motion[33l [33|. In this section we investigate the effects 



C. Effects of the Kerr-down conversion 
nonlinearity on the effective tempreture of the 
mirror 

Here we examine the effects of the Kerr and the gain 
nonlinearities on the back-action ground-state cooling of 
the oscillating mirror. In order to investigate the cooling 
of the mirror it is sufficient to consider the mean number 
of quanta of the vibrational excitation of the mirror as 
defined byfUl 

n m = kBTe ff = ksT ( )3 / 31 s 

hw e ff hhJ m J e ff U)eff 
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FIG. 6. (Color online)The normalized effective frequency 
u e ff/io m versus the ui/uj m (a) for rj = 0.01Hz and for different 
values of the gain parameter: G — (blue dashed curve), and 
G = 10 6 Hz (red solid curve), (b) for G = 
ferent values of anharmonicity parameter 
curve), and rj — 0.01Hz (red solid curve). The parameters are 
Um/2-K = 10 MHz, L = 2mm, m = 10 ng, Q = 5 x 10 5 Hz 
P = 6.9 mW, n = 0.01w m , A = u m , and 6 = tt/2 . 
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FIG. 7. (Color online)The normalized effective damping 
rate of the oscillating mirror r f e ff/'m'y m versus the u)/uj m (a) 
for 77 = 0.01Hz and for different values of gain parameter: 
G — (blue dashed curve), and G = 10 6 Hz (red solid curve), 
(b) for G = 10 6 Hz and for different values of anharmonicity 
parameter: rj — (blue dashed curve), and rj = 0.01Hz (red 
solid curve). The parameters are uj m /2ir — 10 MHz, L = 
2mm, m = 10 ng, Q = 5 x 10 5 Hz, P = 6.9 mW, K = 0.01w m , 
A = uj m , and 9 = 7r/2 . 



where we expand the effective frequency and effective 
damping rate of the oscillating mirror in a Taylor series 
around u> = uj m and keep only the leading terms in the 
respective expansions 

w e//( w ) ~ W e//K) = W e //, 
7e//M ^ 7e//(w m ) = Jeff- (32) 

The ground-state cooling is approached if n m < 1. In 
FigsEJa) and[UJb) we have plotted, respectively, n m and 
T e ff for the initial equilibrium temperature T = 0.4K 
versus A/uj m and for different values of G and rj. It is 
evident that as G increases the system cools down to the 
ground-state (n m < 1) while as rj increases the effective 
temperature of the system increases . The heating of the 
system in the presence of the Kerr medium is due to the 
energy shift of the cavity which weakens the field-mirror 
coupling and reduces the number of intracavity photons 
because of the photon-photon repulsion mechanism. 



V. EFFECTS OF THE KERR-DOWN 
CONVERSION NONLINEARITY ON THE 
INTENSITY AND QUADRATURE NOISE 
SPECTRA OF THE TRANSMITTED FIELD 

The intensity spectrum of the transmitted field is given 
by the Fourier transform of the two time correlation func- 
tions (8al ut (t + T)8a out (t)): 

dT(Sal ut (t + T)6a out (t))e-^. (33) 

-oo 

Solving the matrix equation (| 10[) we obtain the solution 
for the cavity-field fluctuation operator 5a and then using 
the input-output relation a out = \[2na — ai n , we obtain 
the fluctuations of the output field to be 

Sa out (uj) = Vi(w)£(w) + V 2 (uj)Sa in (uj) + V 3 (uj)Sa\ n (uj), 

(34) 
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FIG. 8. (Color online) Plots of (a) the mean number of 
vibrational quanta n m and (b) the effective temperature T e / / 
for G = 2 x 10 7 Hz,77 = 0.05Hz (blue dotted-dashed curve), 
G = 2 x 10 7 Hz,77 = 0.01Hz (red thick curve) and G = 0, r\ = 
0.01Hz (green thin curve). The parameters are u> m /2n = 10 
MHz, L = 10mm, m = 10 ng, Q = 5 x 10 5 Hz, P = 6.9 mW, 
k, = 0.7cj m , and 6 — tt/2 . 



and Sa' out (uj) = [£<w(— , where 



Vi(w) 
V 2 (w) 

V 3 (w) 



= V2k 



■{<2 s (k — iA 



2k ihgLa 2 . 2 



(«i + r)}-i, 

2k ^ff^|a s | 2 2 
i U h i w ™ 



a; 



- iw) - 2Ga*e }, 



d(u) 1 M 
(s + ifw + Ai))}, 



(35) 



In Eq. (|34[) the first term stems from the mechanical os- 
cillator thermal noise, while the other two terms are from 
the input vacuum noise. Making use of the correlation 
peroperties for the noise forces and Eqs. ([M)l and ([55*]) . we 
obtain the output field spectrum as 



S(lu) = \V 3 (-uj)\ 2 + mfvy m u)[l + coth( 



2k B T> 



\Vi(-u)\ 2 . 
(36) 



The quadrature noise spectrum of the transmitted filed 
is given bvj36j 



S»= / dT{5x^{t + T)5x^{r)) ss e-^ 

J — OO 

- (^Hfe^'H), (37) 

where 5x° v ut {uj) = e- lv 5a out (uj) + e^Sal ut (uj) is the 
Fourier trnasform of the output quadrature, with ip as its 
externally controllable phase angle which is experimen- 
tally measurable in a homodyne detection scheme (37| .For 
the system under consideration quadrature noise spec- 
trum is given by 

S V ( U ) = e-**C%?{w) + e 2 ^C a °r M 

+ Cr>)+C11H, (38) 



where 



C%* = {m7l7 TO a;[l + coth( 
+ V 2 (u)V 3 (-u)}, 

C2t={mh lm uj[l + coth( 
+ \V 2 (u J )\ 2 }, 

CTa = {mfi 7m w[l + coth( 
+ \V 3 (-u J )\ 2 }. 



Hid 

2k B T 
2k B T 
2k B T 



JlViH^c-w)) 
)MMI 2 
)]|^(-^| 2 



(39) 
(40) 
(41) 



We define the optimum quadrature squeezing S op t(uj) by 
choosing f{uS) in such a way that dS v {uj)/dtp = 0. This 
yields 



2iip opt _ 



no 



\GZ l \ 



(42) 



Then, substituting back into Eq. (j3"8)) . we obtain the op- 
timized squeezing spectrum as 

S opt {u) = -2\C:f\ + C a 7» + C%!(u>). (43) 

In FigsOHa) and [9jb) we have plotted, respectively, the 
intensity and the optimized squeezing spectra of the 
transmitted field for different values of the gain nonlin- 
earity versus the normalized response frequency u)/u) m . 
As can be seen increasing the gain parameter G causes 
the intensity and the degree of squeezing of the transmit- 
ted field to increase. The shift in the intensity spectrum 
is due to the energy shift of the cavity in the presence of 
the Kerr medium. In FigsfTUJa) and fTOT b) we have plot- 
ted, respectively, the intensity and the optimized squeez- 
ing spectra of the transmitted field for different values of 
the anharmonicity parameter 77 versus the normalized re- 
sponse frequency w/uj m . It can be seen as 77 increases the 
intensity and the degree of squeezing of the transmitted 
field increases. 




FIG. 9. (Color online) (a) The intensity spectrum and (b)the 
squeezing spectrum of the transmitted field versus the normal- 
ized frequency u)/uj m for G = 1 X 10 7 Hz (red thick curve), 
G = 2 x 10 7 Hz(blue dashed curve), and G = 3 x 10 7 Hz( green 
thin curve). The parameters are rj = 0.01Hz, u> m /2n = 10 
MHz, L = 10mm, m = 10 ng, Q = 5 X 10 5 Hz, P = 6.9 mW, 
K — 0.7cj m , and 6 = n/2 . 



FIG. 10. (Color online) (a) The intensity spectrum and 
(b)the squeezing spectrum of the transmitted field versus the 
normalized frequency u)/u m for rj = 0.01 Hz (red thick curve), 
77 = 0.03 Hz(blue dashed curve), and rj = 0.05 Hz(green thin 
curve). The parameters are G = 10 7 Hz, uj m /2-K — 10 MHz, 
L = 10mm, m = 10 ng, Q = 5 x 10 5 Hz, P = 6.9 mW, 
K = 0.7uj m , and 6 — ir/2 . 



VI. ENTANGLEMENT PEROPERTIES OF THE 
SYSTEM 

It has been shown previously [38f that within each 
optomechanical cavity the photon-phonon entanglement 
can be generated by means of radiation pressure . Here 
we are interested in the generation of stationary entan- 
glement between the movable mirror and the cavity field 
in the presence of Kerr-down conversion medium. Since 
our system is in a Gaussian state and the linear nature 
of Eq. lfTO]) preserves the Gaussian nature of the initial 
state of the sy_stem, we can use the logarithmic negativ- 
ity measure [39j to quantify the entanglement. When the 
stability conditions of Eq. (fT^)l are fulfilled, we can solve 
Eq. tfTTjf for the 4x4 stationary correlation correlation 
matrix (CM) V 



MV + VM 1 



D, 



(44) 



where the elements of the correlation matrix V 



are defined as 



UjUi)/2 and D 



Diag[Q,mhrf m uj m {2n + 1),k,k] is the diagonal diffusion 
matrix in which we used the following apporximation 



u) coth( 



tlUJ 







2k B T 



ui m {2n + 1), 



(45) 



2ksT fiw„ 

where n = [ e R " m / fcbT — Then the photon-phonon 

entanglement can be quantified by the logarithmic nega- 
tivity En as 



whe 



re rj 



E N = max[0, -ln2r/"], , (46) 

, r , il/2 

= 2" 1 / 2 T,{V) - v/S(F) 2 - 4detF , is the 

lowest symplectic eigenvalue of the partialtranspose of 
the 4 x 4 CM, V, with £(V) = detV A + detV B - 2detV c , 
and we used the 2x2 block form of the CM 



V A V c 
V b 



(47) 



with Va associated to the oscillating mirror, Vb to the 
cavity mode, and Vc describing the optomechanical cor- 
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relations. In Fig. (|ll|) wc have plotted Em versus the nor- 
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FIG. 11. (Color online) The logarithmic negativity En 
versus the normalized detuning A/uj m for (77 = 0, G — 0) 
(green thick curve), (77 = 0.01Hz,G = 0)(red dashed-dotted 
curve), (17 = 0, G = 6 x 10 6 Hz)(blue dashed curve) and 
(77 = 0.01Hz, G — 6x 10 6 Hz)(orange thin curve). The param- 
eters are ui m /2TY = 10 MHz, L = 1mm, m = 10 ng Q = 5 X 10 5 
Hz, P = 15 mW, k = 0.9o; m , and = tt/2. 
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FIG. 12. (Color online) The logarithmic negativity En ver- 
sus the normalized detuning A/ujm for G = 10 3 Hz(green 
solid curve), G = 10 6 Hz(red dashed-dotted curve), and 
G = 10 7 Hz(blue dashed curve) . The parameters are r\ = 0.01 
Hz, cj m /27r = 10 MHz, L = 1mm, m = 10 ng, Q = 5 x 10 5 
Hz, P = 15 mW, k = 0.9o; m , and = tt/2. 

malized detuning A/uj m for a bare cavity, a cavity that 
contains only the Kerr medium, a cavity that contains 
only the gain medium and a cavity that contains both 
the Kerr and gain medium. We see that, as expected, 
while the Kerr medium reduces the degree of phonon- 
photon entanglement, the gain medium increases it con- 
siderably. This result arises from the fact that the radia- 
tion pressure, which is responsible for the phonon-photon 
entanglement, increases considerably in the presence of 
the gain nonlinearity while decreases in the presence of 
the Kerr nonlinearity. Therefore similarly to what hap- 
pens for the mirror cooling process there is a competition 



between the gain and Kerr nonlinearities for increasing 
and decreasing quantum effects, respectively. In FigJT2l 
we have plotted En versus the normalized detuning for 
different values of the gain nonlinearity. It can be seen 
that as G increases the degree of entanglement increases. 
It should be noted that wc cannot increase G arbitrarily 
because of the stability condition of the system. 



VII. SUMMARY AND CONCLUSIONS 

In this paper we considered the canonical optomechan- 
ical cavity composed of a partiallly transmitting mirror 
and one perfectly reflecting movable mirror containing 
a Kcrr-down conversion nonlinear medium. For this sys- 
tem we investigated the influence of the two nonlinearities 
on the dynamics of the oscillating mirror, the intensity 
and squeezing spectra of the transmitted field and the 
steady-state mirror-field entanglement. For the dynam- 
ics of the oscillating mirror we have found that in the 
peresence of the two nonlinearities, there is a competi- 
tion between the two quantum processes; the splitting of 
the pump photons in two degenerate subharmonic pho- 
tons as well as increasing the intracavity photon number 
due to the nonlinear gain medium, and shifting the en- 
ergy level of the cavity proportional to the anharmonicity 
parameter rj which blocks the entrance of a large num- 
ber of photons into the cavity. This competition can 
lead to enhancement of the radiation pressure contribu- 
tion to the normal modes of the cavity and NMS be- 
low the threshold goa s < k. Also, it was demonstrated 
that in the presence of the Kerr-down conversion nonlin- 
earity the ground-state cooling of the mechanical mirror 
is possible. Although with increasing the anharmonic- 
ity parameter 77 the vibrational excitation and effective 
temperature of the mirror increases, the enhancement of 
the gain medium can compensate the increasing of the 
phonon number to some extent. In the investigation 
of the intensity spectrum and quadrature squeezing of 
the transmitted field, it was shown that the range of re- 
sponse frequency of the field undergoes a shift because 
of anharmonicity parameter 77 but in this range of the 
frequency both the Kerr and gain media increase the de- 
gree of squeezing. Finally, wc investigated the influence 
of the Kerr-down conversion on the degree of station- 
ary entanglement between the cavity and the mechanical 
modes. It was shown that the Kerr nonlinearity reduces 
the degree of entanglement considerably, while the gain 
medium increases it. 
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